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Abstract 

For the scalar Wick-Cutkosky model in the particle representation we perform a similar 
variational calculation for the 2-point function as was done by Feynman for the polaron 
problem. We employ a quadratic nonlocal trial action with a retardation function for 
which several ansatze are used. The variational parameters are determined by mini- 
mizing the variational function and in the most general case the nonlinear variational 
equations are solved numerically. We obtain the residue at the pole, study analytically 
and numerically the instability of the model at larger coupling constants and calculate 
the width of the dressed particle. 
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1 Introduction 



The need of nonperturbative methods in in quantum physics is obvious considering the many 
problems where strong couphng and/or binding effects render perturbation theory inadequate. 
In nonrelativistic many-body physics variational methods are widely used under these circum- 
stances while this is not the case in relativistic field theory. However, Feynman's successful 
treatment of the polaron problem shows that variational methods may also be used for a 
nonrelativistic field theory provided that the fast degrees of freedom can be integrated out and 
their effect properly taken into account in the trial action. In a previous paper (henceforth 
referred to as (I) ) we have extended the polaron variational method to the simplest scalar 
field theory which describes heavy particles ( "nucleons" ) interacting by the exchange of light 
particles ("mesons"). In euclidean space-time the Lagrangian of the Wick-Cutkosky model is 
given by 

C = \ {d,^f + 1mo^$2 + 1 {d,^Y + ImV^ - ^7$V • (1) 

Here Mq is the bare mass of the nucleon, m is the mass of the meson and g is the (dimensionfuU) 
coupling constant of the Yukawa interaction between them. In the quenched approximation 
the meson field can be integrated out and one ends up with an effective action for the nucleons 
only 

/|i,-^e-'-('<--'.'.... (2) 

Note that this is formulated in terms of trajectories x{t) of the heavy particle ("particle rep- 
resentation") which are parametrized by the proper time r and obey the boundary conditions 
a;(0) = and x(/3) = x. To obtain the 2-point function one has to perform the path integral 
over all trajectories and to integrate over (3 from zero to infinity with a certain weight. It is, 
of course, impossible to perform this path integral exactly and, again following Feynman, we 
have approximated it variationally by a retarded quadratic two-time action. In (I) we have 
proposed various parameterizations for the retardation function which enters this trial action 
and derived variational equations for the most general case when its form was left free. 

The purpose of the present paper is to investigate numerically these parametrizations as well 
as to solve the variational equations. This fixes the variational parameters which will be used 
to calculate physical observables in forthcoming applications. One quantity which we evaluate 
in the present paper is the residue on the pole of the propagator. Another one is related to 
the well-known instability of the Wick-Cutkosky model : although the effective action (Q) is 
very similar to the one in the polaron model the ground state of the system is only metastable. 
This does not show up in any order of perturbation theory but, as we have demonstrated in 
(I), the variational approach knows about it. Indeed an approximate solution of the variational 
equations has revealed that there are no real solutions beyond a certain critical coupling. In 
the present paper we will find exact numerical values for this critical coupling and calulate the 
width of the unstable particle for couplings beyond it. 

This paper is organized as follows: The essential points of the polaron variational approach 
are collected in Section |[ while Section ^ is devoted to the numerical methods and results. In 
Section ^ we investigate the instability of the Wick-Cutkosky model in our variational method 
and determine analytically and numerically the width of the dressed particle. The variational 
principle can also be applied away from the pole, which is explored in Section |^ and used to 
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calculate the residue at the nucleon pole. The main results of this work are summarized in the 
last Section. 



2 Polaron Variational Approach 

Following Feynman's treatment of the polaron problem we have performed in (I) a variational 
calculation of the 2-point function with the quadratic trial action 



St[x] 



13 1 

dr -x^ + 
2 



/3 r-ri 2 

dri / dT2f{Ti-T2) [x{ti) - x{t2)] . 
Jo 



(3) 



Here /(ri — T2) is an undetermined 'retardation function' which takes into account the time 
lapse occurring when mesons are emitted and absorbed on the nucleon. In actual calculations 
we rather have used the Fourier space form 



(4) 



k=0 



where the bk are the Fourier components of the path x(r) and the Fourier coefficients Ak are 
considered as variational parameters. The variational treatment is based on the decomposition 
of the action ^eff into ^efr = St + AS and on Jensen's inequality 



(5) 



Near = — Mpj^y^ the 2-point function should behave hke 



G2{P') 



(6) 



where < Z < 1 is the residue. As was shown in (I) this requires the proper time /3 to tend to 
infinity. One then obtains the following inequality 



,,0 M? A 
M\ < ^ + - M\ 



+ - {n + V) 

A 



where 



inA: 



(7) 



(8) 



4:71^ m 

is a finite mass into which the divergence of the self-energy has been absorbed and A a variational 
parameter. For /3 — > 00 all discrete sums over Fourier modes Ak turn into integrals over the 
'profile function' A{E — kn/P) and one finds 



2 r 

TT Jo 



dE 



hiA{E) + 



A{E) 



- 1 



(9) 



as well as 
V = 



^ „ , da du 

STT^ Jo Jo 



/ AMphysO- ^ 
e m//((7), T-^,u 



— — e {m\/a^ 0, u) 
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(10) 
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Here we use the abbreviations 



and 



e(s, t, u) 



I — u 



exp -- 



u 



2" 



4 /"CO 

^/o A{E) 



1 sin2(Ea/2) 



^2 



(12) 



Because /i^(c'") behaves hke cr and a /Aq for small and large a, respectively, we have called it a 
'pseudotime'. Note that in Eq. ( p!0D the particular renormahzation point /^o = has been used 
to regularize the small-cx behaviour of the integrand. As we have shown in (1) the total result 
is, of course, independent of /xq- 

The profile function A[E) is linked to the retardation function /(a) by 



A{E) =i + I da f{a) sin^ . 
In (1) we have studied the following parametrizations 

'Feynman' parametrization: 



(13) 



which leads to 



'Improved' parametrization: 



which entails 



ME) 



fF{a) = 
Af{E) 
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(14) 
(15) 

(16) 
(17) 



In both cases v, w are variational parameters whose values have to be determined by minimizing 
Eq. (0). 

As well as the above parametrizations, it was possible not to impose any specific form for the 
retardation function but to determine it from varying Eq. (|^) with respect to A and A{E). 
This gave the following relations 



vr^ Jo 
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Together with Eq. (JT^ they constitute a system of coupled variational equations which have 
to be solved. Assuming /i^((T) ^ a and m ~ we have found in (I) an approximate solution 
which had the same form of the retardation function as the 'improved' parametrization and 
exhibited the instability of the system beyond a critical coupling constant. In the general case 
we can read off the variational retardation function from Eq. 



/var(c^) 



327r2 ^\a) 



du 



m 



1 + -" 



u 



u 



2/.2(a) 



u 



AM 



phys 



a 



(20) 

Obviously it has the same 1 / cr^-behaviour for small relative times as the 'improved' parametriza- 
tion (|l^). Finally we mention that by means of the variational Eq. (0), one can find the 
following expression for the the 'kinetic term' Q defined in Eq. (^ 



iTC^ Jo 



da 




(21) 



This is demonstrated in the Appendix and will be used in Chapter |[ That the kinetic term 
Q can be combined with the 'potential term' is a consequence of the virial theorem for a 
two-time action |Q which the variational approximation fulfills. 



3 Numerical Results 

In this Section we will compare numerically the various parametrizations for the retarda- 
tion function. Because we are primarily interested in an eventual application in pion-nucleon 
physics, we have chosen the masses and coupling constants appropriately. Of course, the model 
does not really give a realistic description of the pion-nucleon interaction as spin- and isospin 
degrees of freedom as well as chiral symmetry are missing. 
In short, we use 

m = 140 MeV (22) 
Mphys = 939 MeV (23) 

and the results are presented as function of the dimensionless coupling constant 

The relevant quantity for the physical situation is the strength of the Yukawa potential between 
two nucleons due to one-pion exchange 0, which is approximately given by (depending on the 
spin-isospin channel) 



V2 



2 



where g'^/4:'n' = 14 is the pion-nucleon coupling. In the Wick-Cutkosky scalar model the 
corresponding strength is just the dimensionless coupling constant a that we have defined in 
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a/u fMeVl 


^/w MeV] 


A 


Ml fMeVl 


A(0) 


0.1 


1850 


1845 


0.97300 


890.23 


1.0120 


0.2 


1805 


1794 


0.94400 


839.73 


1.0257 


0.3 


1756 


1739 


0.91250 


787.29 


1.0417 


0.4 


1702 


1678 


0.87773 


732.69 


1.0606 


0.5 


1641 


1608 


0.83843 


675.70 


1.0838 


0.6 


1569 


1527 


0.79223 


616.09 


1.1142 


0.7 


1477 


1424 


0.73355 


553.93 


1.1582 


0.8 


1325 


1254 


0.63714 


490.60 


1.2485 



Table 1: Variational calculation for the nucleon self-energy in the Wick-Cutkosky model using 
the 'Feynman' parametrization (|r^) for the profile function. The parameters v,w obtained 
from minimizing Eq. (pTf ) are given as well as A and the intermediate renormalized mass Mi 
(see Eq. (R)). The last column lists ^4(0) = v^/w"^. 



Eq. (^41). It should also be remembered that a Yukawa potential only supports a bound state 
if 



m 
M 



a > 1.680 ^ = 0.2505. (26) 



We have minimized (cf. Eq. (^) 

- Ml < (A^ 



- 2A) M^, + 2in + V) 



(27) 



with the 'Feynman' ansatz (|I5|) and the 'improved' ansatz (|T3). This minimization was per- 
formed numerically with respect to the parameters X,v,w by using the standard CERN pro- 
gram MINUIT. The numerical integrations were done with typically 2 x 72 Gauss- Legendre 
points after mapping the infinite-range integrals to finite range. For the 'improved' retardation 
function we had to calculate /x^((j) and fl numerically. Tables |1] and |^ give the results of these 
calculations. We also include the value of Mi although it doesn't have a physical meaning: 
finite terms (which, for example, arise when a different renormalization point is chosen) can be 
either grouped with Mi or with V. However, from the variational inequality ( pTf ) we see that 
Ml is a measure of the quality of the variational approximation: the larger Mi the better the 
approximation. 

Although the value of the parameters v and w are rather different for the Feynman and the 
'improved' parametrization, the parameter A and the value of the profile function at E = 
are very close. This is also reasonable when we study the behaviour of these quantities under 
a reparametrization of the particle path: it can be shown that a rescaling of the proper time 
j3 ^ 13/k leaves the variational functional invariant if 



K 



We are working in the 'proper time gauge' 
parameters v, w then obviously are different 



(28) 



K = I. 

see Eqs. 



In a general 'gauge' /t the variational 
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^/^; fMevl 


a/u; [Me VI 


A 


A-fi fMeVl 


) 


0.1 


677.2 


674.6 


0.97297 


890.25 


1.0158 


0.2 


661.4 


656.0 


0.94390 


839.78 


1.0338 


0.3 


640.8 


632.3 


0.91223 


787.43 


1.0548 


0.4 


613.7 


601.9 


0.87715 


732.97 


1.0808 


0.5 


596.7 


581.2 


0.83741 


676.20 


1.1109 


0.6 


570.4 


550.7 


0.79040 


616.97 


1.1514 


0.7 


534.3 


509.2 


0.72996 


555.45 


1.2118 


0.8 


468.2 


434.5 


0.62429 


493.44 


1.3482 



Table 2: Same as in Table |l| but using the 'improved' parameterization fH 
function. 



for the profile 



but A{Q) = v'^/w'^ and A are gauge- invariant. 

For both parametrization no minimum of Eq. (|2^) was found beyond 



a > ar 



where 



0.824 
0.817 



'Feynman' 
'improved' 



(30) 



(31) 



This value of the critical coupling is surprisingly close to the value — vr/4 which we obtained 
from the approximate solution of the variational equations in (I). On the other hand when the 
parameter A is fixed to A = 1, i.e. a less general trial action for "momentum averaging" (see 
(I)) is used, then a minimum is found for all values of a. This points to the important role 
played by this parameter. Indeed, in the approximate solution of the variational equations 
found in (I) the branching of the real solutions into complex ones is most clearly seen in the 
approximate solution for A. We can also trace the instability to the inequality (^ for the 
physical mass : a clear minimum as a function of A exists only as long as the coefficient of 
1/A , i.e. Ml/ 2 + VL + V stays positive. However, with increasing coupling Mi shrinks and V 
becomes more negative until at the critical coupling the collapse finally occurs. 

We have also solved the coupled nonlinear variational equations ([T8|), ([T9| ) together with 
(^) numerically Q . This was done by the following iterative method: we first mapped variables 
with infinite range to finite range , e.g. 



E 



a 



1 

taxiijj 



(32) 
(33) 



and then discretized the integrals by the standard Gauss-Legendre integration scheme, with 
typically 72 or 96 gaussian points per integral. The functions A{6), fi'^lip) as given by the 

^Note that the variational sohition is also reparametrization invariant: Eqs. and ([l^ ) are consistent 

with the condition (28). 
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A 


Ml fMeVl 


A(0) 


0.1 


0.97297 


890.25 


1.0151 


0.2 


0.94389 


839.78 


1.0322 


0.3 


0.91223 


787.43 


1.0520 


0.4 


0.87718 


732.97 


1.0755 


0.5 


0.83738 


676.20 


1.1044 


0.6 


0.79030 


616.98 


1.1421 


0.7 


0.72968 


555.47 


1.1972 


0.8 


0.62262 


493.55 


1.3188 



Table 3: The variational parameter A, the renormalized mass Mi and the value of the profile 
function at E = from the solution of the variational equations . 

variational equations were then tabulated at the gaussian points using as input the values of 
X, A{6), ^'^{ip) from the previous iteration. We started with the perturbative values 

aw = A{e,p = i 



f^'i^if'^ = ]^^tan7/>, (34) 



1 

phys 

and monitored the convergence with the help of the largest relative deviation 



A„-Max|^ AW ' A{e:jM ' Ji^^^) j' ^"l'^- 

(35) 

Some numerical results are given in Table ^. Comparing with Table |^ we observe a remarkable 
agreement with the values from the 'improved' parametrization. It is only near a = 0.8 that 
the variational solution is appreciably better as demonstrated by the numerical value of Mi 
which measures the quality of the corresponding approximation. 

This may also be seen in Figs. |I|, |^ and ^ where the different profile functions and 
pseudotimes are plotted for a = 0.2 and a = 0.8. One can also confirm from the graphs that 
the numerical results indeed have the limits for a, E either small or large which we expect 
from the analytical analysis. Furthermore, it is clear that the 'improved' parametrization of 
the trial action is in general extremely close to the 'variational' one, while the 'Feynman' 
parametrization deviates much more. Finally, it is interesting to note that the profile function 
of the 'variational' calculation has a small dip near E = which is a result of the additional 
terms in the retardation function (|20|). These rather innocent looking deviations will become 



important if an analytic continuation back to Minkowski space is performed in which physical 
scattering processes take place. 

Examples for the convergence of the iterative scheme are shown in Fig. ^ . It is seen that 
for small coupling constant we have rapid convergence which becomes slower and slower as the 
critical value 

ttc = 0.815 ( 'variational' ) (36) 



8 



1.035 
1.030 
1.025 
^ 1.020 

Ld 

^ 1.015 

1.010 
1.005 
1.000 



1 1 

- \ 


1 1 1 1 1 
a = 0.2 


\ 


\ 




- ■ ■ ■ - \ 




\ • 




y ■ . . 




1 1 


1 1 1 1 1 







3 



4 5 
E/M^ 



Figure 1: Profile function A{E) as function of E for the 'Feynman' parameterization (|T5D 
(dotted line), the 'improved' parameterization (|1^ (dashed line) and the 'variational' solution 
(solid fine). The dimensionless coupling constant is a = 0.2. 



1.000 




2 4 6 8 

Figure 2: Ratio of pseudotime /i^(cr) to proper time a for a = 0.2 . The labeling of the curves 
is as in Fig. |l] . An expanded view of the small-cr region is shown in the inset. 
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Figure 4: Convergence of the iterative solution of the variational equations as a function of the 
number of iterations n. The convergence measure A„ is defined in Eq. (^). 
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1.3 




100 200 300 400 500 

m (MeV) 

Figure 5: Critical coupling constant as a function of the meson mass m. The nucleon mass is 
fixed at M = 939 MeV . The crosses indicate the points at which the critical coupling has been 
calculated, the line through them being drawn to guide the eye. 

is reached. Finally, beyond a > only a minimal relative accuracy can be reached and the 
deviations increase again with additional iterations. 

How the critical coupling depends on the meson mass is shown in Fig. |^ . It turns out 
that the good agreement of the approximate value of ~ vr/4 with the numerical value 
obtained for m = 140 MeV was an accidental one: at m = we have = 0.641. There is 
also a surprisingly strong but nearly linear m-dependence which we cannot reproduce from 
an approximate solution of the variational equations when taking m ^ but still assuming 
/i^(cr) ^ a. 

4 Instability and Width of the Dressed Particle 

In all parametrizations of the profile function A{E) which we investigated numerically in the 
previous section it turned out to be impossible to find a (real) solution of the variational 
equations or the variational inequality for coupling constants above a critical value ac- This is 
a signal of the instability of the model which is already seen in the classical "potential" 

V^'\^, = ^Mo^$2 + ImV^ - 9^' ^ (37) 
and tells us that the physical mass of the dressed particle becomes complex 

Mphy. = M - i^. (38) 

In the following we take the real part of the nucleon mass to be M = 939 MeV and try to 
determine the width F. 

Note that in a perturbative calculation no sign of the instability shows up : the one-loop 
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self-energy 

9^ 

is perfectly well-behaved. Also the one-loop effective potential 0, |^ is not very indicative: in 
quenched approximation it is given by 

A detailed analysis shows that the quantum corrections lower the barrier which makes the 
ground state metastable in V^^\^, (p) but do not remove it. In addition, the one-loop effective 
potential develops an imaginary part but it is easy to see that Im v}^^ vanishes for — 2(7$ > 
Ag'^cfP' /m?, i.e. it contains a (non-analytic) step function. Therefore all proper one-loop Green 
functions (which are generated by the effective action) carry no sign of the instability. 
In contrast the variational approach for the two-point function knows about the instability 
if we allow the parameter A in the trial action to vary. Since the approximate solution of 
the variational equation for A in (I) clearly showed the impossibility to obtain a real solution 
beyond ac , we will first study the width of the state using similar approximative methods 
before turning to the exact numerical evaluation. 



1 + 



■u + 



Mi 



u 



m? 1 



u 



(39) 



m 



2 p2 + Mi - 2gip - ie 



(40) 



4.1 Approximate analytical treatment 

In order to discuss complex solutions of the variational equations it is useful to introduce the 
complex quantity 

C = AMphys (41) 

and to write it in the form 

C = Co e-'^ . (42) 

It is a phase x 7^ which will lead to the complex pole (|38D of the two-point function. With 
the same approximation m = and ^^{c) ~ a which was used before we now obtain 

^ ' (43) 



A 



a M2 



Here the dimensionless coupling constant is defined in terms of the real part of the physical 
mass 



9^ 



47r2M2 ■ 

Due to Eq. (plD the kinetic term vanishes under the same approximation 



(44) 
(45) 



and the potential term becomes 

^2 ^1 ^oo \ 
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du 
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exp 
m2 \ — u 
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exp 



m 1 — u 



-a 



u 



(46) 



12 



These are rather drastic simphfications but the exact numerical calculations show that the 
imaginary part of Q is indeed smaller (by a factor of five) than Im V. Note that V is not infrared 
stable, i.e. it diverges if the meson mass m is set to zero. With the above approximations the 
stationarity equation (0) then reads 



1 + 



u 



m? 1 — u 



du In 



1 + 



1 — u 



(47) 



(4^ 



Using Eq. (^3]) this is equivalent to 

o o 2a; 

(2 = A/r + 

TT IT 

If we take the imaginary part of this equation it is possible to set m = and we obtain 

Co' sin2x= — M^x. 

How do we determine the width of the unstable state ? We take the defining equation (|4l|) for 
( , ehminate A by means of Eq. (H3[) and use Eq. (pSf). This gives 



(49) 



M 



r ^ aM^ 
'2 = « + — 



(50) 



The real and imaginary parts of this equation allow us to express Co and the width as a function 
of the phase x- A simple calculation gives 



Co 



M 



2 cosx 



1 



'1 



4a 



vr 



cos^ X 



and the width is 



2Mtanx 



1 4a 2 
1 cos^ X ■ 

TC 



(51) 



(52) 



We have chosen the root which results in a positive width for < x < 7r/2. Finally, substituting 
Eq. ( |5TD into Eq. ( ^Qj ) gives the transcendental equation which determines the phase x- After 
some algebraic transformations we obtain it in the form 



a 



TT 



2x sin2x 



(2x + sin 2x)2 cos^ x 



(53) 



It is easy to see that the function h{x) grows monotonically from h{0) = 1 to h[Ti/2) = oo. 
Solutions x(q^) therefore only exist for 



a > Or 



TT 
4 



(54) 



which is the same critical value of the coupling constant at which previously the real (approxi- 
mate) solutions of the variational equations ceased to exist. It is also easy to find solutions for 
the transcendental equation (ESf) for small x • from h{x) = 1 + x' + C^(x^) '^6 find 



X 



la — Ur 



Olr 



(55) 
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where, of course, Oc = 7r/4 should be used. Since the expression (|^) for the width can be 
transformed into 

T = 2Mtanx- — , . „ 56 
we obtain the following nonanalytic dependence of the width on the coupling constant 

M -) . (57) 



3 \ Oc / 
This should be valid near the critical coupling constant 



4.2 Numerical results 



For the numerical solution of the complex variational equations we follow the approximate 
analytical solution as closely as possible. However, some of the relations used previously do 
not hold exactly. For example, the quantity 



(^2 ^oo g-2 / 

L = — / do-——— / du u e \ nifiicr) 
2 JO (o") JO \ 



u 



(5^ 



would be unity for m = 0, /U^(cr) = a but has some complex value in the exact treatment. Sim- 
ilarly, 7^ and V deviate from the aproximate value (|46| ) . Without invoking the simplifying 
assumptions Eq. (H) changes to 



(2 = Ml - —M^L + 2{n + V) 
Following the same steps as in the approximate treatment we obtain 

M 



(59) 



Co 



which replaces Eq. ( pi]) and 



2 cosx 



a = Tc 



4a 

1 + Wl cos X Re {Le'^) 



K 

Re {Le'^) + K cosx 



(60) 



(61) 



which supersedes Eq. (|53|). Here 



K 



sin2x 



Im 



(62) 



Instead of Eq. ( ^6]) one can show that the width itself has now the exact form 

K sinx — Im (Le*'^) 



2M 



K COSX + Re {Le^^) 



(63) 
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Table 4: The width F of the unstable state from the complex solution of the variational 
equations for a > = 0.815. The width is given as a function of the phase x which 
determines the corresponding coupling constant a according to Eq. (^Tj). The complex value 
of the profile function at = is also listed. 

We have solved the coupled complex equations by specifying a value for the phase x and 
determining the corresponding value of the coupling constant a by means of Eq. ([61|). Of 
course, this could be done only iteratively by starting with 

= 1, 7^(0) ^2x/ sin 2x, 
= A(o)(E) = l. 

Typically 20 - 25 iterations were needed to get a relative accuracy of better than 10~^. Table ^ 
gives the results of our calculations. It is seen that the width grows rapidly after the coupling 
constant exceeds the critical value. In Fig. ^ this is shown together with the approximate 
(small-x) behaviour predicted by Eq. (0). After the critical coupling constant in this formula 
has been shifted to the precise value one observes a satisfactory agreement with the exact 
result. 

Finally Figs. |^ and ^ depict the complex profile function A{E) and the complex pseudotime 
/i^((T) for X = 0.5 , i.e. a = 1.153 . Compared to the real solutions below Oc ( cf. Figs. |l] - H ) 
one does not notice any qualitative changes in the real part of A{E) as one crosses the critical 
coupling. 

5 The Two-point Function Away from the Pole 

Up to now we only have determined the variational parameters on the nucleon pole. However, 
the variational principle also applies to ^ ~^phvs- This forces us to consider sub-asymptotic 
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Figure 6: Width of the unstable state as a function of the couphng constant as obtained from 
the solution of the complex variational equations (see Table ^ ). The dashed line shows the 
approximate solution 
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Figure 7: Real part of the profile function and of the ratio of pseudotime to proper time for 
a = 1.153. 
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Figure 8: Imaginary part of the profile function and of the ratio of pseudotime to proper time 
for a — 1.153. 



values of the proper time (3. We first deal with the residue at the pole which gives us the 
probability to find the bare nucleon in the dressed particle. 



5.1 The residue 

To calculate the residue it is most convenient to use the "momentum averaging" scheme de- 
veloped in (I) because in this approach there are only a few subasymptotic terms. To be more 
specific the quantity p,'^{a,T) introduced in Eq. (1.98) has an additional term which exactly 
cancels the l/f3 term which arises from application of the Poisson summation formula. With 
exponential accuracy we therefore have 



4 roo 

fi\a,T) ~ - / dE 

TT Jo 



1 sm^{Ea/2) 



A{E) 



E^ 



(64) 



This is a big advantage as we do not have to expand the potential term <^ 5*1 S> in Eq. (1.97) 
in inverse powers of (3. The only source of subasymptotic terms in <^ 5*1 ^ is then from 
the T- integration from a/2 to (3 — a/2 which simply gives a factor (3 — a. Applying the 
Poisson formula to the kinetic term Q defined in Eq. (1.100) we obtain, again with exponential 
accuracy 



_ 9 roo 

n(p) ^ - dE 

IT Jo 



In A{E) + 



A{E) 



- 1 



+ 



InA(O) + 



A{0) 



- 1 



We recall from (1) that the 2-point function may be written near the pole as 

1 f°° 

G2{p) - 7^ / dp exp 

2 Jo 



(65) 



(66) 



Collecting all non-exponential terms the function F[j3,p'^) therefore has the large-/? expansion 



F(/5,p2) ^ Fo(p^) + ^i^i(p'), 



(67) 
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where 



p2 + M^-p^{i-\Y + n 



3_ 

47r^ Jo 



da j duel m^{a), 



—iXpa 



u 



(6^ 



is what we have used before on the nucleon pole {p = iMphys) and 



Fi{p^) = \nA{0) + 



1 - ^(0 ) , 9 
A{0) 



+ 



^7r^ Jo 



da 



a 



fi^[a) Jo 



du e (mn(a), — , u 



(69) 



Note that the potential term in -Fo(p^) develops a small-a singularity which renormalizes the 
bare mass Mq but -Fi(p^) is finite. 

Neglecting the exponentially suppressed terms and performing the proper time integration 
we thus obtain the following expression for the two-point function 



exp 



lnFo(j9')-Fi(p^ 



It is now very easy to calculate the residue Z at the pole (see Eq. 
the point p^ 



(70) 

by expanding around 



-Mpjjyg where Fq vanishes. We obtain 



exp 



-Fi(-M2 



physJ 



^o(-MpV) 

where the prime denotes differentiation with respect to p^ . Explicitly we find 

iVoiVi 



D 
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(72) 



No 
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exp -lnA(O) + 1 
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u 



(73) 
(74) 

(75) 



In the last line the stationarity Eq. (18) for A was used to simplify the denominator D. Note 
that this also applies to the case where one parametrizes the profile function A{E). This 
demonstrates that 

Z = ^ (76) 
A 

is always positive. It seems to be more difficult to prove in general that Z < 1 although all 
numerical calculations clearly give this result. Finally, it is again useful to check the variational 
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Ct 


' Ppvti man' 


'irnnrovpd ' 


'variational' 


nPTtnrl)ativp 

l-J V_jX U l-XX V-J CAj U X V V_j 


0.1 


0.96090 


0.96087 


0.96087 


0.96200 


0.2 


0.91934 


0.91914 


0.91918 


0.92399 


0.3 


0.87467 


0.87418 


0.87428 


0.88599 


0.4 


0.82600 


0.82494 


0.82521 


0.84798 


0.5 


0.77184 


0.76996 


0.77036 


0.80998 


0.6 


0.70940 


0.70610 


0.70672 


0.77198 


0.7 


0.63216 


0.62597 


0.62697 


0.73397 


0.8 


0.51086 


0.49187 


0.49284 


0.69597 



Table 5: Residue at the pole of the two-point function for the different parametrizations of 
the profile function. The heading 'Feynman' gives the result in the Feynman parametrization 
whereas 'improved' refers to the improved parametrization from Eq. (p^. The residue calcu- 
lated with the solution of the variational equations is denoted by 'variational'. For comparison 
the perturbative result is also given. 

result in perturbation theory. With A{0) = 1 + 0{g'^) one sees that Nq = 1 + 0{g^). Similarly 
(1 — A)^ = 1 + 0{g'^). Expanding A'^i and 1/A to order g"^ we obtain 

Z = 1 [ da [ du (1 — u) exp f 

87r2 io io ^ ' ^ \ 2 u 

87r2 io Mp\y,M2 + m2(l - + "^^^ ' " 

This coincides with what one obtains from the perturbative result for the self-energy (^) in 
the usual way. 

Table ^ contains the numerical values of the residue obtained with the different parametriza- 
tions as well as the perturbative result from Eq. (|77D 

^perturb = 1 - 0.38004 « . (78) 

It is seen that for a near the critical value appreciable deviations from the perturbative result 
occur. For example, at a = 0.8 perturbation theory says that there is a probability of nearly 
70% to find the bare particle in the dressed one whereas the variational results estimate this 
probability to be less than 50%. It should be also noted that the residue is not sua infrared 
stable quantity, i.e. for m — > Z also vanishes. From the variational equations one can deduce 
that 

Z const, (79) 

with K = a/(7rA^). For massless mesons the residue at the nucleon pole must vanish because it 
is well known (e.g. from Quantum Electrodynamics) that in this case the two-point function 
does not develop a pole but rather a branchpoint at p"^ = — Mpi^ . 



-^.)+0(.^) 

(77) 
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5.2 Variational equations for the ofF-mass-shell case 



It is also possible to apply the variational principle away from the pole of the two-point function 
by varying Eq. (|70D . This gives 



^0) 



Note that on mass-shell where Fq vanishes the previous variational equations follow. We will 
not elaborate on Eq. (|^) further but only point out that the perturbative self-energy (pOj) is not 
obtained from the off-shell variational equations (PP]). In the limit fi'^{cr) a, A{E) — > 1, A — > 1 
one rather finds 
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47r2 



, A' 9 
In — + 



47r2 



du In 



u 



1 — M 



9' 



47r2 



duu- 



p^ + Ml 



u 



^1) 



which is an expansion of Eq. (^) around p^ = — Mq (or Mf which is the same in lowest 
order perturbation theory). The reason for this somehow unexpected result is the neglect 
of exponentially suppressed terms in deriving Eq. ([7^). Indeed it is easy to see that one 
obtains the correct perturbative self-energy only if the upper limit of the a-integral is kept 
at (3 and not extended to infinity as we have done in deriving Eq. ([70|) . The difference is 
one of the many exponentially suppressed terms which we have neglected. Thus, the off- 
shell variational equations (|80D only hold in the vicinity of the nucleon pole and in order 
to investigate variationally the two-point function far away (say near the meson production 
threshold = — (Mphys+?7i)^ ) one has to include consistently all terms which are exponentially 
suppressed in j3. This is beyond the scope of the present work. 



6 Discussion and Summary 

In the present work we have performed variational calculations for the 'Wick-Cutkosky polaron' 
following the approach which was developed previously [0. We have determined different 
parametrizations as well as the full variational solution for the retardation function which 
enters the trial action. Since the nucleon mass is fixed on the pole of the 2-point function the 
value of the functional which we minimize is of no physical significance but only a measure of 
the quality of the corresponding ansatz. This is in contrast to the familiar quantum- mechanical 
case where an upper limit to the ground-state energy of the system is obtained. However, our 
calculation fixes the variational parameters with which we then can calculate other observables 
of physical interest. 

One of these quantities was the residue on the pole of the propagator for which we have 
compared numerically the results of the variational calculations to first order perturbation 
theory in Table table: residue. For small couplings all results for the residue agree, since in 
this case the variational approach necessarily reduces to perturbation theory independent of the 
value of the variational parameters. What is rather remarkable is that for larger couplings the 
three parametrizations of the profile function in our variational approach yield rather similar 
results, which are now of course different from the perturbative calculation. As we have seen, 
the 'improved' and 'variational' actions have the same singularity behaviour, for small relative 
times, as the true action, so here one might expect some similarity in the results. This is 
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however not true for the 'Feynman' parameterization which has a rather different form, so its 
agreement with the other two is not preordained. This similarity is also exhibited in Tables |l| 
and for A and ^(0), but of course not for the parameters v and w which enter the respective 
profile functions and which are 'gauge' (i.e. reparametrization)-dependent quantities. Also 
the critical coupling at which real solutions ceased to exist was nearly identical in all three 
parametrizations. The similarity of the results for the different ansatze presumably indicates 
that these results are not too far away from the exact ones. 

We were not only able to determine the critical coupling but also to deduce qualitatively 
and quantitatively the width which the particle acquires beyond the critical coupling. This was 
achieved by finding complex solutions of the variational equations first approximately by an 
analytic approach and then exactly by an iterative method which closely followed the analytic 
procedure. Although the present approach does not describe tunneling (which we expect to 
render the system unstable even at small coupling constants but with exponentially small width 
) the polaron variational method is clearly superior to any perturbative treatment in this 
respect. 

We have concentrated mostly, although not exclusively, on the on-shell 2-point function, 
i.e. the nucleon propagator. This corresponds to the limit where the proper time goes towards 
infinity. It is possible, however, to go beyond the on-shell limit. This was necessary, for example. 



for the calculation of the residue of the 2-point function in Section 5.1. Nevertheless, the residue 



is a quantity which is calculated at the pole and thus only requires off-shell information from 
an infinitesimal region around it. This has the effect that the variational parameters for the 
calculation of the residue are the same as the on-shell ones. As one moves a finite distance away 
from the pole the variational parameters themselves become a function of the off-shellness 
(see Section |5^ ). 

In conclusion, we think that the present variational approach has yielded nonperturbative 
numerical results which look very reasonable and are encouraging. We therefore believe it 
worthwhile to try to extend it in several ways. First, in a sequel to this work we will generalize 
the present approach to the case with n external mesons and thereby study physical processes 
like meson production or meson scattering from a nucleon. This can be done by employing the 
quadratic trial function whose parameters have been determined in the present work on the 
pole of the 2-point function. Such a 'zeroth order' calculation is similar in spirit to a quantum 
mechanical calculation in which wave functions determined from minimizing the energy func- 
tional are used to evaluate other observables. More demanding is the consistent 'first-order' 
variational calculation of higher-order Green functions as this requires the amputation of pre- 
cisely the non-perturbative nucleon propagators which have been determined in the present 
work. That this is indeed possible will be demonstrated in another paper in this series. 

Of course, finally we would like to apply these non-perturbative techniques to theories which 
are of a more physical nature. Among these one may mention scalar QED, the Walecka model 



10, 11 1 or QED. The latter two will require introduction of Grassmann variables in order to deal 



with spin in a path integral. As such, this should not pose a fundamental problem. A greater 
challenge, however, is to extend such an approach beyond the quenched approximation or to 
nonabelian theories where the light degrees of freedom cannot be integrated out analytically. 
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Note added 

After completion of this work we became aware of the pioneering work by K. Mano |12[ in which 
similar methods are applied to the Wick-Cutkosky model with zero meson mass. Mano uses the 
proper time formulation, the quenched approximation and the Feynman parametrization for 
the retardation function to derive a variational function for the self-energy of a scalar nucleon 
(the expression following his Eq. (6.18)) which is identical with our Eq. (0) after proper 
identification of quantities is made. However, for minimizing the variational function Mano 
sets (in our nomenclature) v = w{l + e), expands to second order in e and finds an instability 
of the ground state for glig^^^/8nM'^ > 0.34. Note that guano = 9 so that this translates 
into a critical coupling ~ 0.22 which is much smaller than the value which we obtain from 
the exact minimization. In addition, in the present work we consider non-zero meson masses, 
employ more general retardation functions, and calculate residue and width of the dressed 
particle. 
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Appendix : An alternative expression for Q 



var 



Here we derive Eq. (^T]) for the kinetic term Q when the variational equations are fulfilled. 
We first perform an integration by parts in the definition of Q. The slow fall-off of the 
variational profile function with E 
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We then write the variational equation ([T9|) for A{E) in the form 
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The integration over E can now be performed giving a factor 7r/i^(cr)/4 due to Eq. (p!^). 
Therefore we have 
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which is just one term in the expression (^3) for VL. To get the other one we differentiate the 
variational equation for A{E) with respect to E and observe that 
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One has to be careful not to interchange the i?-integration and the a-differentiation. We 
therefore perform an integration by parts and obtain 
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Note that the boundary term at a = gives a contribution because of X(0) = 1. This 
contribution exactly cancels the term (yf^/Svr^ in Eq. ([A.2|) . Combining both terms for Q 
(which do not exist separately due to the slow fall-off of A{E) ) we obtain 



dE 



j^ A'jE) ^ 1 - AiE) 
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from which Eq. (0) follows. In the last line again an integration by parts has been performed 
but this time there is no contribution from the boundary terms. 
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